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ABSTRACT 

The  relative  distance  between  various  Antisubmarine  Warfare  vehicles 
is  an  important  parameter  for  the  evaluation  of  sensor  equipment.   This 
thesis  presents  a  computer  model,  utilizing  a  time  dependent  error  navi- 
gation system,  Omega,  for  the  determination  of  this  relative  distance 
parameter.   A  methodology  for  predetermining  a  relative  distance  error 
distribution  for  a  specific  area  by  simulating  various  Antisubmarine 
Warfare  vehicle  tracks  from  known  navigational  errors  has  been  estab- 
lished.  Attempts  to  smooth  estimated  tracks  of  Antisubmarine  Warfare 
vehicles  which  receive  the  same  time  dependent  navigational  information 
simultaneously  was  found  not  to  be  advisable  as  this  disturbs  the  correla- 
tion of  navigation  error  between  the  vehicles. 
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Section  One 
Introduction  to  the  Problem 
1-1  General  Discussion  of  the  ASW  Exercise 

The  fleet's  Antisubmarine  Exercises  are  an  excellent  example  of  the 
operational  environment  in  which  test  and  evaluation  of  weapon  systems  can 
be  conducted.   These  exercises  are  necessary,  for  the  sophistication  and 
cost  of  the  present  state-of-the-art  of  military  weapon  systems  requires 
that  their  capability  be  determined  in  this  test  and  evaluation  phase. 

The  ASW  exercises  not  only  provide  the  opportunity  to  observe  the 
new  system  but  also  to  train  personnel  and  determine  the  effectiveness  of 
the  existing  tactics  and  equipment.   However,  this  method  of  evaluation 
is  expensive  which  places  additional  emphasis  on  the  essentiality  of  ob- 
taining the  best  results  which  will  in  fact  reflect  the  true  capability 
of  the  weapon  system. 
1-2  Problem 

One  of  the  primary  purposes  of  an  ASW  exercise  is  to  provide  data  to 
evaluate  the  potential  of  detecting  devices  of  the  ASW  forces.  The  inter- 
play of  the  different  vehicles  involved  (i.e.  aircraft,  surface  ships  and 
submarines)  in  an  ASW  exercise  create  a  unique  problem  in  that  not  only 
are  their  exact  positions  never  established  but  also  confirmation  of  de- 
tection by  actual  sightings  are  not  always  possible.   This  situation 
breeds  the  continually  arising  question,  "Were  any  two  vehicles  within  de- 
tection range  of  the  detecting  sensors?" 

This  question  is  hopefully  resolved  during  the  formal  reconstruction 
of  the  ASW  exercise  which  takes  place  after  its  termination.  The  person- 
nel usually  represented  at  these  sessions  are  a  mixture  of  those  directly 
involved,  along  with  the  necessary  analyst  and  administrative  personnel. 


One  of  their  first  tasks  is  to  determine  the  tracks  of  the  various  vehicles 
by  reconstructing  them  from  the  respective  navigators'  logs.   The  present 
reconstruction  procedure  is  spelled  out  in  great  length  which  results  in 
a  process  that  is  tedious  and  time  consuming.  |_lj  Unfortunately,  even  after 
these  tracks  have  been  reconstructed,  serious  doubt  exist  to  their  accu- 
racy.  The  cause  of  this  doubt  does  not  arise  from  the  reconstruction  pro- 
cedures but  basically  in  its  source  of  navigational  information — the 
navigator's  log. 

The  navigator's  log  is  composed  of  information  gathered  from  available 
navigational  equipment  on  board  the  respective  vehicles.   The  navigational 
errors  present  in  the  log  are  the  result  of  an  accumulation  of  navigation- 
al equipment  and  human  errors.   The  relative  distance  between  any  two  ve- 
hicles in  question  is  of  primary  interest  in  this  reconstruction  session. 
The  error  in  this  relative  distance  between  any  two  vehicles  (i.e.  see 
figure  1)  results  from  a  combination  of  their  separate  inherent  naviga- 
tional errors.   Knowledge  concerning  this  relative  distance  error  is 
presently  not  well  formulated.   But  this  relative  distance  error  is  one 
of  the  parameters  necessary  in  analyzing  the  effectiveness  of  the  detec- 
tion equipment. 
1-3  Objective 

General  Electric  TEMPO,  Santa  Barbara,  California  is  presently  under 
contract  with  the  United  States  Navy  to  investigate  the  problems  intro- 
duced by  the  relative  distance  errors  and  to  propose  means  of  alleviating 
or  reducing  them.  Smoothing  techniques  are  being  considered  which  would 
attempt  to  reduce  the  existing  relative  errors  by  smoothing  the  estimated 
positions  found  from  present  navigation  information.  The  employment  of  a 
new  type  of  navigational  aid  called  OMEGA,  which  is  in  the  operational- 


Antenna^B"  (observed  position) 
Antenna »-Bw  (  known  position  ) 


°^S  0900 


estimated  track  for  Destroyer- 
actual  track  for  Destroyer. 


Anter_na"A"  (observed  position  ) 

G   S02C0 
Antenna"AI!  (known  ' 


QQ9Q0 


relative  distance  of  actual 
positions  at  time  090 0, 

•ive  distance  of Est- 
imated positions  ai 
time  0900. 


Eelative  Distance  Error —Relative  distance  of  actual  position; 

.   ative  distance  of  estimated  position; 


-Rel- 


Omega  receiver  error  =  Known  positions  of  antenna's  -Observed  pos- 
itions for  antenna's. 


Figure  # 


4 


evaluation  stage,  is  under  consideration  as  the  primary  source  of  naviga- 
tional information. 

It  is  therefore  our  objective  to  investigate  the  aforementioned  alter- 
natives and  recommend  a  procedure  that  will  aid  considerably  in  the  track 
reconstruction  problem.   Our  recommendations  will  focus  on  two  specific 
aspects,  first  a  faster  and  more  effective  method  of  reconstructing  the 
various  ASW  vehicles  tracks  and  second  to  provide  a  means  by  which  one 
could  conduct  a  statistical  error  analysis  of  the  tracks. 

Hopefully  our  solution  would  effectively  free  the  personnel  at  the 
reconstruction  session  from  doing  the  actual  track  reconstructing  them- 
selves, thus  permitting  them  to  devote  more  time  and  energy  towards  strict- 
ly analyzing  the  exercises. 
1-4  Approach 

The  first  step  is  to  tackle  the  basic  underlying  problem,  that  of 
errors  introduced  by  all  the  various  methods  in  which  the  navigation  in- 
formation is  presently  obtained.   It  cannot  be  overemphasized  that  if  er- 
rors must  exist,  it  is  highly  desirable  that  these  errors  be  dependent  in 
respect  to  the  navigational  systems  of  both  vehicles,  when  relative  accu- 
racy is  of  interest.   It  follows  that  one  should  then  concentrate  on  find- 
ing the  optimal  combination  of  navigational  equipment  required  to  provide 
this  relative  distance  information.   Clearly  the  minimum  navigational  in- 
formation necessary  is  the  position  of  each  vehicle  at  a  given  time.   How- 
ever, in  order  that  analysis  be  conducted,  a  sequence  of  such  points  must 
exist  so  that  a  continual  picuture  of  the  interplay  is  displayed.   The 
time  interval  between  the  sequence  of  points  depends  on  the  speed  at  which 
the  vehicles  are  traveling.   The  time  interval  between  aircraft  may  be 
one  minute  while  for  submarine  and  surface  vessels  every  five  minutes  may 
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be  sufficient  depending  on  the  detail  desired  of  the  reconstructed  tracks. 
This  sequence  of  points  is  then  basically  a  "dead-reckoning"  track  between 
navigational   fixes   for  each  vehicle. 

The  present  methods  employed  to  obtain  the  DR  tracks   are  to   begin  at 
a  fix  obtained  from  a  navigation  source  and  apply  course  and  speed  infor- 
mation as  one's  guide  until   another  fix  is  established.     This  requires 
continuous  monitoring  of  compasses   and  speed   indicators   if  one  expects  to 
obtain  the  required  accuracy  desired  for  detection  equipment   evaluation. 
This  process  has  not  proven  desirable  even  after  highly  sophisticated 
equipment  has  been  employed  to   automatically  plot  the  tracks.      One  may 
suggest  that  more  fixes  be  obtained   so  that  the  DR  time  can  be   shortened 
thus  resulting  in  better  accuracy.     At  present  this  suggestion  is  sound  in 
theory  but   is  not  practical.      The  reason  being  that   short  time   intervals 
(e.g.    less  than  fifteen  minutes)   between  fixes   for  a  majority  of  present 
Naval    crafts    is  not  easily  obtained. 

The   idea  of  receiving  continuous  fixes   is   intriguing,    so   let  us  pur- 
sue  it   further.      A  position  can  be  obtained   in  two  ways,    first  by  knowing 
the  range  and  bearing   from  a  reference  point.      The  second  method   is  to 
find  two   lines  of  positions  which  intersect.      Since  reference  positions 
(i.e.    land  masses)   are  non-existant   the  majority  of  the  time,   the   inter- 
section of  two    lines  of  position  appears  to  be  the  more  desirable  tech- 
nique. 

With  this   in  mind,    if  one  could  monitor  on  board  the  ASW  vehicles, 
some  sort  of  navigational    input   signals  that  would  establish  their  posi- 
tions, then  a  continuous  fix  could  be  realized.     This  could  be  accomplished 
by  putting  on  board  the  vehicles  recorders  which  record  these  signals  and 
the  corresponding  time.      This  procedure  would  also  take  the  place  of  the 
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navigator's  log  -  eliminating  human  errors. 

Hopefully,  if  errors  are  introduced  using  this  technique,  they  will 
be  dependent  to  this  technique  and  will  not  affect  appreciably  the  relative 
accuracy  between  any  two  vehicles. 

The  optimal  navigation  system  selected  to  be  used  in  conjunction  with 
this  technique  is  one  that  is  not  only  compatible  to  a  continuous  readout 
but  also  possessing  errors  which  result  from  a  source  common  to  all  ve- 
hicles in  the  immediate  vicinity,  e.g.  variations  in  the  velocity  of  wave 
propagation.   We  propose  that  the  Omega  navigation  system  now  in  the  op- 
erational-evaluation stage,  which  possesses  these  desired  characteristics, 
be  selected.   However,  any  navigational  system  with  these  characteristics 
may  be  used. 

In  fact,  we  now  have  a  navigation  system,  consisting  of  Omega  navi- 
gational information,  which  can  be  monitored  continuously  and  simultan- 
eously by  all  vehicles. 
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Section  Two 
Omega 
2-1  Introduction 

Omega  is  a  Very  Low  Frequency  (VLF)  navigation  system  operating  in 
the  internationally  allocated  navigation  band  in  the  electromagnetic  spec- 
trum between  10  and  14  Kc/s.   (See  Table  1  for  further  system  character- 
istics) .  |_2J  As  in  Loran  each  pair  of  stations  provides  the  navigator  with 
a  family  of  hyperbolic  lines-of-position.   The  5000  to  6000  nautical  mile 
baseline  between  stations  gives  a  global  coverage  capability  with  a  maxi- 
mum of  only  eight  stations. 
2-2  Propagation  Medium 

Most  modern  radio  navigation  systems  are  designed  to  obtain  optimum 
results  from  propagation  medium  involved.   From  this  follows  the  fact  that 
most  navigation  systems  are  propagation  limited.   Omega  is  no  exception. 
Capabilities  and  operational  procedures  stem  in  large  part  from  the  prop- 
erties of  VLF  propagation.   It  should  be  noted  that  the  propagation  prop- 
erties to  be  considered  are  solely  those  at  the  Omega  frequencies.   Much 
data  has  been  acquired  at  the  basic  10.2  Kc/s  frequency  so  definite  con- 
clusions can  be  made.   Omega  is  a  simple  system  with  range  and  accuracy 
based  on  the  features  of  VLF  propagation,  second  order  effects  are  well 
understood  although  large.   Easy  to  use  corrections  are  prepared  for  use 
with  Omega. 

The  most  apparent  feature  of  an  Omega  VLF  phase  difference  measurement 
as  observed  at  some  remote  stationary  point  is  a  diurnal  variation.  While 
day-night  transitions  result  in  large  changes  in  position,  they  are  nor- 
mally quite  regular  and  predictable.   If  the  diurnal  variation  can  be  pre- 
dicted it  need  only  be  subtracted  from  observed  data  to  yield  a  navigation 
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Receiver  Specification    (Ground  Monitor) 

Frequency:      10.2   KC   fixed  or  continuously 
tunable   from  8.2   to    15.2   KC . 

Sensitivity:     0.75  db  noise  figure. 

Differential   Amplitude:      80  db. 

Input   Impedance:      50  ohms. 

Detection:      Synchronous. 

Instrumental  Accuracy:      0.2%  phase-difference 

of  two   noise-free   signals. 

Tracking:      Zero  velocity  error   servos. 

Signal   Conditions:      Track  reliably   in  signal-to- 
noise  ratio   of   1/10. 

Commutator  Synchronization:      Semi-automatic   operation 

in  signal-to-noise  ratios 
to    1/3. 

Reference  Oscillator:      Long-time   stability   1   part   in 

10°  per  day,    short   time  stability 
5  parts    in   10  "  per    30  minutes. 

Size:      Height  57  %",  Width   19  %",    Depth   17  %'*. 

Weight:      486  pounds. 

Power:      115  volts,    single  phase,    60   cycles, 
12   amperes. 


Table   1 
Omega  Receiver  Specifications 
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variate  which,  to  a  good  approximation,  varies  quite  regularly  with  posi- 
tion. 
2-3  Accuracy 

Factors  effecting  the  accuracy  (i.e.  the  magnitude  of  the  navigation- 
al error)  to  be  obtained  from  a  practical  navigation  system  include:  re- 
peatability of  observations  such  as  day-to-day  variations  at  a  particular 
site,  accuracy  with  which  observed  readings  may  be  associated  with  real 
terrain  features,  and  degradation  experienced  in  actual  operations  due  to 
shipboard  installation  or  procedures. 

Although  Omega  is  an  absolute  system,  relative  accuracy  is  important 
in  many  applications.   Because  the  limitations  on  Omega  accuracy  is  the 
phase  stability  of  a  transmission  through  long  distances,  phase  perturb- 
ations of  a  signal  from  a  remote  transmitter  at  any  two  receiving  sites 
in  the  same  neighborhood  will  be  much  the  same.   Relative  accuracy  is 
therefore  significantly  better  than  absolute  accuracy.   It  is  also  indica- 
tive of  relative  positioning  capability  with  a  task  force.   As  separation 
distance  between  receivers  decreases,  the  relative  accuracy  increases. 

The  additional  difference  between  the  experimental  and  operational 
systems  must  be  taken  into  account.   The  experimental  system  used  short- 
term  synchronization,  that  is,  propagational  fluctuations  on  the  path  be- 
tween two  Omega  transmitters  which  were  fed  into  the  phase  radiated  by  one 
of  the  transmitters.   In  the  operational  system,  the  synchronization  will 
be  achieved  through  the  use  of  very  stable  oscillators,  using  a  time  con- 
stant of  one  or  more  days.   Thus  the  phase  will  be  radiated  at  a  constant 
value  and  corrected  only  slowly  to  keep  it  at  the  assigned  value. 

The  effect  of  this  procedure  would  reduce  the  operational  error  to 
2/3  those  of  the  experimental  system.   Thus  we  can  predict  that  the 
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standard  deviation  of  the  operational    line-of-position  errors   for  day- 
time will  be  0.32  miles   and  nighttime,   0.82  miles.    [_3J 

In  the  case  of  the  no n- correlated   lines  of  position,   the  root-mean- 
square  fix  error   is   usually  taken  as  the   square  root  of  the   sum  of  the 
squares  of  the  two    line-of-position  errors   divided   by  the  sine  of   the 
crossing  angle  between  the  two    lines  of   position.      If   the  two   position 
lines  have  the   same  accuracy  we  can  compute  a  daytime  rms   fix  error  of 
0.48  miles   and   a  nighttime  rms   fix  error  of   1.22  miles.    j~3j 
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Section  Three 
Determination  of  Relative  Distance  Error  Distribution 
3-1  Omega  Implementation 

At  the  initial  planning  stages  of  an  ASW  exercise,  recorders  capable 
of  recording  Omega  information  from  the  receivers  which  in  turn  receive 
signals  simultaneously  from  two  different  pairs  of  pre-selected  Omega 
transmitters  will  be  installed.   The  tapes  of  raw  recorded  Omega  data  for 
each  individual  vehicle  will  be  transferred  to  a  computer  facility  when 
convenient  or  at  the  termination  of  the  ASW  exercise.   The  tapes  will  then 
be  processed  in  order  that  the  navigation  information  is  in  proper  form 
for  reading  into  a  computer  program.   A  computer  program,  such  as  the  one 
referenced  in  Appendix  B,  will  reconstruct  the  vehicle  tracks  in  such  a 
manner  as  to  provide  the  user  with  relative  distance  information  of 
selected  vehicles  (i.e.  permutations  of  aircraft  and  destroyers  with  sub- 
marines). With  this  relative  distance  information,  an  analysis  regarding 
whether  or  not  any  two  vehicles  were  within  detection  range  for  a  partic- 
ular sensor  can  be  accomplished  with  the  aid  of  a  relative  distance  error 
distribution. 
3-2  Preliminary  Determination  of  a  Relative  Distance  Error  Distribution 

The  actual  tracks  along  with  the  estimated  tracks  of  the  vehicles 
must  be  available  in  order  to  determine  the  relative  distance  error  dis- 
tribution for  an  area.   The  actual  track  of  a  specific  vehicle  is  defined 
as  that  track  which  is  not  in  error.   The  estimated  track  is  defined  as 
that  track  which  incorporates  the  error  inherent  to  the  Omega  navigation 
system.   With  the  above  two  tracks  defined,  the  relative  distance  error 
at  any  given  time  is  that  value  obtained  by  subtracting  the  distance  be- 
tween the  actual  positions  of  the  two  vehicles  from  the  distance  between 
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their  estimated  positions,  as  shown  in  figure  1. 

It  is  immediately  apparent  that  the  actual  tracks  of  the  various  ve- 
hicles are  not  obtainable  because  a  navigator  cannot  be  sure  of  his  exact 
position.   Our  solution  to  this  aforementioned  problem  has  been  to  simu- 
late hypothetical  ASW  tracks  and  refer  to  them  as  the  actual  tracks  of  the 
ASW  vehicles.   The  estimated  tracks,  on  the  other  hand,  corresponding  to 
these  actual  tracks  can  be  determined  by  incorporating  known  errors  from 
a  pair  of  Omega  receivers.   These  Omega  receiver  errors  can  be  obtained 
by  observing  the  difference  between  the  known  fixed  position  of  the  re- 
ceiver's antenna  and  the  position  observed  from  the  resulting  navigational 
information  being  received  at  a  given  time,  as  depicted  in  figure  1.   It 
follows  that  these  actual  tracks  will  be  constructed  in  the  general  area, 
a  reasonable  distance  from  the  antenna,  and  in  a  corresponding  time  frame 
with  the  antenna.   At  this  point  it  is  necessary  for  our  determination  of 
the  relative  distance  error  distribution  to  assume  that  the  dominate  error 
inherent  to  the  Omega  navigational  system  is  in  fact  the  variation  in  the 
velocity  of  the  electromagnetic  waves  in  the  media  and  is  not  appreciably 
affected  by  the  mode  of  travel  of  the  receivers.   This  assumption  is  based 
on  the  fact  that  the  Omega  system  operates  on  the  phase  difference  con- 
cept which  is  not  disturbed  by  the  mode  of  travel,  although  the  signal-to- 
noise  ratio  may  be  affected. 
3-3  Simulation  inputs 

The  simulated  estimated  tracks  make  use  of  Omega  strip  charts  re- 
ceived from  the  Naval  Electronics  Laboratories,  San  Diego,  during  the  per- 
iod of  1-31  March  1963  (see  figure  2).   Strip  charts  were  used  from  three 
different  receiving  antennas  located  in  the  San  Diego  area  and  two  more 
receiving  antennas  located  in  the  Pt.  Mugu  and  Pt.  Huenume  areas.   This 
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gives  a  capability  of  getting  simultaneous  Omega  navigational  information 
from  the  five  different  sources  to  that  simulations  of  various  vehicles 
can  be  achieved. 

To  determine  the  Relative  Distance  Error  Distribution,  a  computer 
program  was  written  to  gather  this  necessary  information  using  as  its 
primary  inputs  raw  Omega  data,  time  and  Longitude- Latitude  data  of  the 
hypothetical  tracks.   It  was  necessary  to  physically  interpolate  the  raw 
Omega  data  from  the  strip  charts.   It  may  be  pointed  out  here  that  this 
interpolation  technique  may  introduce  errors  independent  of  the  Omega  nav- 
igation system  error  itself,  but  the  maximum  error  was  100  yards  in  our 
interpolation  procedure.   It  can  be  assumed  that  if  our  Relative  Distance 
Errors  are  acceptable,  when  a  more  sophisticated  technique  is  employed, 
better  Relative  Distance  Errors  should  result  because  of  the  exclusion  of 
these  independent  errors. 
3-4  Determination  of  Antenna  Errors 

Omega  strip  chart  information  for  two  fixed  antenna  positions  as  a 
function  of  time  is  used  as  input  to  the  computer  program.   The  minimum 
interval  of  time  used  for  the  Omega  readings  was  one  minute  because  of  the 
"seaman's  eye"  technique  used  for  acquiring  this  information  from  the  strip 
charts.   Diurnal  variation  corrections  were  then  applied  to  the  raw  Omega 
data.   This  was  accomplished  by  reading  into  the  program  a  matrix  correc- 
tion table  which  corresponded  to  the  applicable  period  of  time  and  loca- 
tion (i.e.  month,  day,  hour,  and  longitude-latitude).   Since  this  matrix 
varied  not  only  in  time  but  also  in  longitude-latitude,  it  was  necessary 
to  employ  a  linear  interpolation  technique  to  determine  the  appropriate 
correction  for  each  input  reading  (see  figure  3).   As  each  reading  was 
corrected  it  then  went  through  a  longitude- latitude  conversion  routine 
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SKYWAVE  LINEAR  INTERPOLATION  CORRECTION  ROUTINE 

Selected  Latitude-Longitude  Outlay 
50N-40W  50N-36W  50N-32W 

8 — (Vehicle's   last  position) 
46N-40W  46N-36W  46N-32W 

42N-40W  42N-36W  42N-32W 

Omega   station  pair  correction  values   for  Latitude-Longitude 
which  enclose  Vehicle's   position  for   specified  time 

Time-0900 
50  51 

t       I 

49-< 49.5 ^50.5 


1 


B  --(Vehicle's   last 

position)  5^j 


Linear   interpolation  between  values 


Time  Interpolation 

Time  0900-48 
Specified  Latitude-Longitude  value 

Time  1000-53 

Time  of  position  desired  -  0930 

Time(0900-0930)X(55-48)Value  =    3.5 

POSITION  CORRECTION  ~  49.5  +  3.5  =  53 


Figure  3 

Skywave  Linear   Interpolation  Correction 
Rout  ine 
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which  is  outlined  in  detail  in  Appendix  C-III.   Finally  the  vector  error 
was  found  by  comparing  this  converted  reading  which  represents  the  estima- 
ted position  of  the  antenna  with  the  actual  antenna  position  so  that  the 
inherent  error  associated  for  that  particular  antenna  can  be  determined. 
This  process  was  repeated  for  the  entire  input  readings  for  both  antennas. 

At  this  point  it  is  possible  to  determine  the  relative  error  necessary 
to  determine  the  Relative  Distance  Error  Distribution  of  the  reconstructed 
tracks  if  one  is  completely  convinced  that  maneuvering  of  the  ASW  vehicles 
have  no  effects  on  the  relative  error  and  if  possible  smoothing  techniques 
were  not  to  be  investigated.   This  relative  error  can  be  simply  found  by 
algebraically  subtracting  the  relative  distance  of  the  two  actual  positions 
of  the  antennas  from  the  relative  distance  of  their  estimated  positions  at 
a  given  time  (see  figure  1). 
3-5  Hypothetical  Track  Simulation 

To  add  intuitive  appeal  to  the  problem,  actual  tracks  simulating 
typical  ASW  tracks  for  aircraft,  destroyers,  or  submarines  were  developed. 
This  also  allows  the  opportunity  to  investigate  the  effects  of  smoothing 
techniques  on  the  overall  relative  distance  error  distribution.   An  aux- 
iliary program  (i.e.  see  Appendix  A)  was  written  to  accept  Longitude-Lat- 
itude information  for  constructing  an  "actual  track".   To  these  longitudes 
and  latitudes  at  a  given  time  were  added  the  components  of  the  vector 
error  from  one  of  the  antenna's  to  obtain  the  estimated  tracks  (i.e.  for 
example,  antenna  A's  errors  construct  a  submarine's  estimated  track  and 
antenna  B's  errors  construct  a  destroyer's  estimated  track  as  shown  in 
figure  1.) 

Smoothing  techniques  were  then  employed  in  an  effort  to  smooth  the 
derived  estimated  tracks.   First  by  the  use  of  the  technique  of  fitting 
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the  estimated  track  by  the  method  of  Least  Squares  Curve  fitting  with  or- 
thogonal polynomials  (refer  to  Appendix  C-II) ,  and  secondly  by  smoothing 
smaller  segments  of  the  estimated  tracks  known  as  the  technique  of  Expo- 
nential smoothing  (refer  to  Appendix  C-I).  A  new  relative  Distance  Error 
was  determined  by  comparing  this  smoothed  estimated  track  with  the  hypo- 
thetical track  as  a  function  of  time. 

The  overall  computer  program  not  only  has  the  capability  of  printing 
out  all  of  the  pertinent  information  such  as  the  Relative  Distance  Errors, 
time,  longitude  and  latitude  as  shown  in  Appendix  A,  but  also  prints  a 
graphical  representation  of  the  actual,  estimated  and  smoothed-estimated 
tracks  as  shown  in  Figure  4. 
3-6  Results 

Sample  points  used  for  determining  the  Relative  Distance  Error  Dis- 
tribution were  generated  from  the  raw  Omega  data.   Partial  sampling  of  the 
data  was  accomplished  by  using  permutations  of  date  and  time  intervals  from 
the  available  data.  As  an  example,  a  day's  transition  periods  (i.e.  sun- 
rise and  sunset)  were  included  as  well  as  continuous  sections  of  day  and 
nighttime  readings.   The  antennas  (i.e.  those  antennas  located  in  San  Diego, 
Pt.  Mugu,  Pt.  Huenume  area),  which  provided  the  raw  Omega  data,  were  loca- 
ted so  that  trials  could  be  conducted  on  pairs  of  antennas  which  were  not 
only  in  close  proximity  but  also  of  considerable  distance  from  one  another 
(i.e.  Pt.  Huenume  and  San  Diego). 

A  cumulative  Distribution  of  the  Relative  Distance  error  (defined  to 
be  the  probability  that  an  observed  value  of  the  Relative  Distance  Error 
is  less  than  or  equal  to  some  value)  was  found  from  using  the  data  which 
included  all  the  trials.   The  distribution  is  found  to  be  Normal  with  a 
mean  of  -120  yards  as  shown  in  Figure  5.   This  mean  implies  that  the 


23 


cc 
cc 
o 


CO 

cr 


H 


o 

O 


en 

i 


ta  oo  ^ 

73  -C)  / 

I-i  '-<     s \ 

a  <a  K s 

o  o  > — ^ 

O  O  I  ,     I 

H  ii  ' — L- 


\ 

\ 

(J3 

cO 

i — 

cO 

1— 

£3 

>--* 

TZ" 

3  3 


o 
cc 
cr 


<o    6?    v ; 


io   CI 

d  aJ  z: 


Pi  euro  //  4 


cr 


24 


J — l_ 


i  Total  '  (JTPfflUa. 


- I 


ative  Distribution  of  ;A11 


Antennae 


i  I 


i±ii 


m 


' 


. 


■ 


m 


:    , 


• 


. 


__ 


! 


M- 


;  i 


LLL 


I 


• 





-LL 


!    ! 


< 


■ 


I 


M 


. 


H 


_ 


— 


:■    ffl 


_ 


33 


m 


■ 


U-IJ 


■ 


• 


: 


' 


Tit 


• 


■  '  i 


L 

tta 


■  i 


; 

: 


B 


i 


I 


• 


. 


I      h 


■ 


IE, 


■ 


■■ 


■    : 


. 


| 


' 


.         ■ 


: 


1 r 





■ 


d 





± 


• 


• 


~~ 


i 


— 


. 





i 


i : 


_ 


- 


1 


i 


i 


p 


- 


. 


& 


■ 


. 


i 


__ 


VA-U 


: 





' 


f- 


LL 


5 


m 


■ 


.  L 


J_T 


3rjrs; 
r 


__ 


- 


t3 


. 


Li  .| 

is 


m 


4+ 


M 


fl 


-4 


Tf 


:<! 


4-44 


—h 


i  th 


_l  — 


hfn 


: 


t±_t 


H±l 


: 


. 


f 


Figure  ,i;5 

Standard  deviation  r  499.2 

Mean  z  -120  yards 


±E3 


S 


-H-H 


4- 


m: 


ffi 


-H-t 


fi 


_i 


m± 


4-H- 


,_u 


;  '  I 


••: 


__ 


— - 


-4— 


!   : 


! 


J_L 


Tf 


l  i 


i 


I 


*& 


US. 


«£ 


MIL 


$11^4 


i^5 


L.ii2 


iiLtk. 


i  i 


! 


I 


mi 

■?■■■{  •*■■'+ 


w : 


-U-l 


LLii 


;   I 


•ffi-jyiffi 


"5  :     EE     : 


lailes- 


tffl 


II 


41 


; 


i  ■ 


relative  distance  for  the  estimated  tracks  was  on  the  average  120  yards 
longer  than  the  relative  distance  of  the  actual  tracks. 

Investigated  next  was  the  cumulative  Relative  Distance  Error  Distri- 
bution for  the  relative  error  sample  points  associated  with  those  antennas 
located  only  in  the  San  Diego  area.   The  analysis  of  this  distribution 
showed  that  normality  could  not  be  assumed,  as  shown  in  Figure  6,  and  the 
mean  was  3.2  yards.   The  above  two  distributions  (i.e.  Figures  5  and  6) 
cover  the  minimum  and  maximum  distance  for  any  two  vehicles. 

When  the  smoothing  techniques  were  used  on  the  estimated  tracks,  the 
Relative  Distance  Errors  were  about  ten  times  as  large  as  those  found  by 
using  no  smoothing.   These  results  were  attributed  to  the  high  degree  of 
time  dependency  of  the  errors  associated  with  the  estimated  track  and  the 
inability  of  the  smoothing  techniques  to  follow  precisely  the  vehicle's 
maneuvers.   The  comparative  results  obtained  are  shown  in  Figure  7.   It 
is  recommended  that  no  smoothing  features  be  incorporated  when  two  ve- 
hicles are  obtaining  simultaneous  readings  from  their  respective  Omega 
receivers. 
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Section  Four 
Other  Omega  Applications 
4-1      Introduction 

If  the  Omega  navigational    system  were  employed  as  the  primary  data 
source  for  the  reconstruction  of  ASW  vehicles'    tracks,  one  might  expect 
instances   in  which  continuous  track   information  could  not  be  obtained. 
This  could  be  due  to  mechanical  difficulties  with  the  equipment,   or,    in 
the   case  of  submarines  engaged   in  evasive  maneuvers,   running  at  depths  be- 
low the  maximum  depths   for  Omega  reception.      This  requires    investigation  of 
techniques,    applied  to  Omega  data,   to    improve  absolute  navigational   accur- 
acies. 
4-2     Omega-No nOmega  Relative  Error  Accuracy 

when  only  one  of  the  two  vehicles   can  receive  Omega  information  con- 
tinuously,  the  track  of  at   least   this  one  vehicle  should  be  constructed, 
so  that  the  tracks   latitude  and  longitude  are  available  for   short   incre- 
ments of  time.      A  representative  sample  of  the  absolute  errors   associated 
with  these  latitudes   and   longitudes   are  shown   in   Figure  8,    under  the  col- 
umn for  non-weighted  vehicles.      It   should  be  emphasized  that  these  abso- 
lute errors  were  obtained  from  data  collected   in  1963  and  present  day 
refinements  of  Omega  have  predicted  that   standard  deviations  of  the  opera- 
tional  line-of-position  errors  during  the  daytime  will  be  0.32  miles  and 
nighttime,   0.82  miles.    [J3j  Trials   simulating  this   situation  were  conducted 
using  vector  errors  to  determine  the  estimated  track,   which  may  be  depend- 
ent with  respect  to  other  errors   for  that  specific  track,    but   independent 
of  the  errors   associated  with  a  second  estimated  track.      One  of  the  esti- 
mated tracks  was   assumed  to  be  determined  by  the  use  of  Omega  and  the  other 
estimated   track  by  some  other  navigational   source   such  as   satellite,   celes- 
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trial,  loran  or  dead  reckoning.  Omega- No nOmega  trials  were  used  to  deter- 
mine if  smoothing  the  Omega  estimated  track,  can  improve  the  relative  dis- 
tance errors.   This  would  be  expected  if  the  errors  in  the  estimated 
position  for  each  vehicle  are  independent.   Figure  9  shows  the  results 
from  such  trials  and  as  expected  it  appeared  that  smoothing  was  in  fact 
an  improvement.   It  is  suggested  that  the  exponential  smoothing  method  be 
used  under  Omega-No nOmega  circumstances  because  of  the  fact  that  there  are 
no  navigational  error  correlations  between  the  vehicles.   Exponential 
smoothing  was  favored  over  least  squares  because  of  its  ability  to  follow 
the  vehicles  maneuvers  was  more  precise  when  a  high  smoothing  constant, 
say  CX  =  .04,  was  used.   In  actual  practice  the  second  track  would  be  re- 
constructed in  such  a  manner  as  to  obtain  the  best  absolute  accuracy  pos- 
sible. This  may  be  accomplished  by  various  methods  such  as  the  procedures 
now  being  used  in  the  formal  reconstruction  sessions  or  from  recording 
continuous  navigational  information  from  a  secondary  source,  e.g.  ships 
inertial  navigation  system. 

4-3  Improvement  of  Absolute  Accuracies  by  Determining  an  Expected  Esti- 
mated Position. 

Improvement  in  the  relative  distance  error  can  be  realized  under 
these  Omega-No nOmega  circumstances  if  one  could  reduce  the  absolute  error 
of  the  Omega  navigational  system.   Appropriate  weighting  of  other  naviga- 
tional information  is  proposed.   An  expected  estimated  position  is  found 
by  using  other  navigational  information  to  a  particular  ship  in  conjunction 
with  Omega.   This  may  be  accomplished  by  weighting  the  obtained  latitudes 
and  longitudes  by  some  factor,  say  A.,  where  i  denotes  a  particular  type 
of  navigational  system,  and  the  sum  of  the  A.'s  must  equal  one.  The  value 
of  A.  is  obtained  by  a  comparative  judgement  of  the  respective  accuracies 
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longitude 
of  the  navigational    sources    (i.e.      .3   (   iat£tU(je     ^   °^  0mega  plus 

loran 
j    (  longitude  ^  determined  by  celestrial  ).   From  this  expected  estimated 
latltude  satellite 

position  and  the  estimated  position  found  when  using  the  Omega  navigation- 
al system,  a  vector  difference  can  be  obtained.  This  vector  difference  is 
then  applied  to  a  number  of  past  and  future  Omega  positions.   A  sample  of 
the  results  obtained  when  using  this  technique  is  displayed  in  Figure  8 
with  o(  -   0.4.   It  may  be  desirable  to  incorporate  a  time  factor,  say  B. 
to  the  vector  difference  constant  in  order  to  take  into  account  a  de- 
creasing effect  on  the  vector  error.   It  is  recommended  that  this  weight- 
ing technique  be  considered  when  Omega-No nOmega  situations  arise. 
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Section  Five 
Operational  Program 

The  knowledge  gained  from  Sections  Three  and  Four  was  incorporated 
into  the  operational  program  (i.e.  the  program  used  at  the  reconstruction 
sessions)  as  shown  in  Appendix  B.   In  Figure  10  a  flow  diagram  is  refer- 
enced which  displays  the  basic  features  and  the  structure  of  the  opera- 
tional program. 

The  implementation  of  this  operational  program  during  the  formal  re- 
construction of  ASW  exercises  would  provide  relative  distance  information 
when  Omega  is  recorded  continuously  by  both  vehicles.   In  conjunction  with 
this  relative  distance  information,  employment  of  a  previously  determined 
Relative  Distance  Error  Distribution  (i.e.  such  as  the  ones  mentioned  in 
Section  Three)  can  be  used. 

The  implementation  of  this  Relative  Distance  Error  Distribution  will 
be  illustrated  through  an  example  shown  in  Figure  11.   It  is  when  the 
vehicles  come  within  detection  range  of  a  particular  sensor  that  is  of 
primary  concern.   Therefore  of  interest  are  those  sample  points  which  lie 
in  one  tail  of  the  Relative  Distance  Error  Distribution,  which  reflects 
the  situation  when  the  actual  relative  distance  is  less  than  the  estimated 
relative  distance.   For  this  example  the  Relative  Distance  Error  Distribu- 
tion found  for  the  antennas  located  only  in  the  San  Diego  area  will  be 
used,  see  Figure  6. 

The  mean  of  this  Relative  Distance  Error  Distribution  is  so  small 
(i.e.  3.2  yards),  it  is  assumed  for  all  practical  purposes  to  be  zero. 
As  shown  in  Figure  11  positions  (a)  and  (b)  reflect  the  estimated  posi- 
tions of  vehicles  A  and  B  respectively  at  time  0900.   The  Relative  Dis- 
tance Error  Distribution  could  be  superimposed  over  one  of  these  positions, 
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say  (b),  with  the  mean  situated  at  (b).  The  distribution  then  represents 
the  sample  points  of  the  various  estimated  relative  distance  readings. 

It  may  be  desirable  to  know  what  the  relative  distance  between  the 
ships  must  be  so  that,  say  ten  percent,  of  the  Relative  Distance  Error 
sample  points  fall  outside  of  the  assumed  maximum  detection  range  of  a 
particular  sensor.  For  this  example,  relative  distance  between  the  ships 
is  found  to  be  -  468  yards  from  the  maximum  detection  range  as  shown  in 
Figure  6.   Therefore  as  shown  in  Figure  11  if  the  maximum  detection  range 
was  1500  yards,  then  the  estimated  relative  distance  was  10  32  yards. 
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Section  Six 
Summary 

It  was  shown  that  the  position  coordinates  for  a  particular  vehicle 
can  be  obtained  as  a  function  of  time  and  that  a  continuous  track  of  the 
vehicles  may  be  graphically  displayed.   This  information  requires  less 
than  five  minutes  of  CDC  1604  computer  time  for  16  hours  of  continuous 
interplay. 

A  methodology  for  predetermining  a  relative  distance  error  distribu- 
tion was  demonstrated  and  an  actual  distribution  for  the  San  Diego  area 
was  established.   This  was  accomplished  by  using  Omega  navigation  informa- 
tion collected  from  receiving  antennas  whose  exact  positions  are  known. 

Attempts  to  smooth  the  vehicles  estimated  tracks  did  not  improve  the 
relative  distance  errors.   This  apparently  resulted  from  disturbing  the 
navigational  error  correlation  between  the  vehicles.   However,  when  abso- 
lute accuracy  of  navigational  positions  was  desired  (i.e.  when  only  one 
vehicle  was  receiving  Omega  information)  smoothing  the  estimated  tracks 
was  an  effective  feature.   Improvement  on  the  absolute  accuracy  was  also 
observed  when  an  expected  estimated  position  was  determined  by  using  other 
navigational  sources  in  conjunction  with  Omega. 

It  is  recommended  that  in  view  of  these  desirable  characteristics, 
that  the  United  States  Navy  continue  the  investigation  into  the  use  of 
Omega  triangulat ion  information  for  track  reconstruction  techniques  in 
ASW  reconstruction  sessions. 
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269.9650  ' 

269.96667 

269.97033 

269.9770*'' 

269.98867 

269.95833 

269.95500 

269.95867 

269.96533 

269.9650^ 

269.97667 

269.98833 

269.98000 

269.96667 

269.95633 

269.95700 

269.97667 

269.98333 

269.9620  , 

269.96867 

269.91333 

269.96600 

269.96667 

269.96833 

269.97000 


102.08:  :• 
102.07283 

102.08967 

102.079  50 

102.084  3  3 

102.069  17 

102.C6400 

102.06885 

102.06867 

102.06850 

102.068  5  3 

102.08317 

102.07800 

102.087  8  3 

102.0  77  67 

102.C5750 

102.0  873  3 

102.07717 

102.07200 

102.07183 

102.07667 

102.076  50 

102.066  3  3 

102.07617 

102.07600 

102.06583 

102 .0  65  67 

102.0  75  50 

102 .065  3  5 

102 .06517 

102.07000 

102.07583 

102».  06  667 

102.068  50 

102.076  3  3 

102.06417 

102.07600 

102  .('6533 

102.06367 

102.06350 

102.07533 

102.06317 

102.06500 

102.07483 

102.06267 

102.0  64  50 

102.07933 

102.06717 

102.072  00 

102.06383 

102.06667 

102.06350 

102.07133 

102.07017 

102.07800 

102.0  7..  8  3 

102.07067 

102.06550 

102.06" 3  3 

102.0  6rM7 

102.07000 


64.96  59  5 
67.25242 
62.61072 

72.86  85  7 

77.371 10 

72.61955 
73.79231 

73.10  381 
68.94338 
65.26473 
65.76965 
6  1.33788 
6  3.6  5  32:') 
60. 5734 i 
57.76717 
64.55  289 
57. 1928C 
63.45005 
58.83602 
60.65441 
50.12554 
49. 16  573 
58.70073 
57.  .1650 
61.49217 
65.15846 
65.65  052 
65.06106 
59.^.!  1835 
59.46  938 
64.81040 

5  9.98  4  42 
62.  73-; 62 
61.30366 
57.81  118 
63.99473 
61.49213 
65.75  5  67 
66.75164 
68.59023 

6  8.5S742 
63.77751 
62.58  9  34 
60.  ;7235 
65.71769 
65.10  733 
63.69012 
70.78  326 
67.01907 
65.75567 
62.21799 
63.34234 
66.26560 
65.75  168 
59.97  108 
64. 18718 
50.90792 
65.07949 
66.70747 
67. 13541 
64.81040 


.  32803 

.03067 

. 32646 

.02772 

.025  8  3 

.030  8  3 

.'"'32  1  S 

.05C06 

.03160 

.03243 

.03237 

.02881 

.029  8  7 

.027  52 

.03161 

.03601 

.•J2860 

.  0  30  1  8 

.03310 

.032  61 

.03182 

.".'3547 

.03500 

.03235 

.03105 

.03528 

.03321 

.03033 

.03522 

•03512 

.03208 

.03217 

.03365 

.:"'3343 

.03203 

.0  34:- 9 

.03105 
.03375 
.03357 
.03322 

.029  6  6 
.03446 
.05427 
.03179 
.034  12 
.•05370 
.029^5 


117.25   " 

1 17.25307 
1 17.25338 
1 17.24937 
1 17.24635 
1 17.25042 
11 7.250 19 
1 17.2  5018 
1 1 7.2  5256 
1 17.2  5478 
1  17.2  54  49 
1 17.25503 
1 1 7.2  544  7 
1 17.2  547  3 
1 17.2  5307 
1 17.25668 
1 17.25670 
1 17.25470 
M 7.25831 
1 17.25719 
1 17.2  580  1 
1 1 7.2  644- 
1 17.25939 
117. 2  5882 
1 17.2  5603 
! 17.25519 
1 17.2  5490 
1 17.25399 
! 17.2  5934 
1 17.2  5905 
1 17.25486 
1 17.25730 
1 17.2  566  3 
1  17.25728 
1  17.:.5827 
1 17.25615 
1 17.25605 
1 17.25508 
1 17.2  544  6 
1 17.25535 
1 17.25204 
1 17.25643 
1 17.25709 
1 17.25702 
1 17.25524 
1 17.25o39 
1 17.2  542  6 


.03172  1 17.25164 
.2  30  97  1 17.2  553 
.03375  1 17.25508 

I  J. '.  "2   O  /"    1  *1  "7    r\  r   *  **    * 


033  8u 
.054  4  7 
.  0  515  4 
.-290  5 
.034 
.03197 
.0  36  72 
.0  33  40 
.034  5  9 
.03453 
.03208 


1  17.25696 
1 17.25667 
1 17.2  5582 
1 17.25314 
1 17.25664 
t 17.25513 
1 17.2  6436 
1 17.25523 
1  1  7.2.-54  8  9 
1 17.2  5460 
1 17.25486 


32.73017 
32.73306 

32.72  827 
3  2.73126 
3.2.72998 
3  2.73420 
32.73567 
32.734  30 
32.73426 
3  2.73423 
32.73429 

35.73  ...5 
52.73154 
52.72874 
32.73151 
32.73729 
3  2.728  81 
32.73177 
32.73311 
32.73320 
52.75179 
32.73161 
32.73463 
3  2.73191 
32.73205 
32. 73 i     7 
3  2  .  7  .ti  5  0  3 
32.73227 
32.734  96 
32.73502 
32.73380 
32.73262 
32.73468 
32.73414 
3  2.73188 
32.73541 


10 


05 


32.  7„ 

32.73554 

32.73561 

3  2.735  70 

3  2.752  3  9 

32.73565 

32.73514 

32.73235 

5  2.73  5  87 

32.735  54 

52.751 17 

32.73472 

32.73329 

32.73554 

32.73467 

32.73558 

32.75346 

52.73126 

32.73146 

32.73355 

52.7352  7 

32.73506 

32.73654 

32.73660 

32.73380 
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TRACK    SIMULATION    OUTPUT    DATA 


VEHICLE    ONE 


H7. 98735 
117.98231 
117.97577 
117.97275 

117.97409 

117.96174 

117. 95442 

117.9470  9 

117.93837 

117.93032 

117.92361 

117.91284 

117.90558 

117.90014 

117,89264 

117.88526 

117.87855 

117.87284 

117.86679 

117.85991 

117.85180 

117.84303 

117.83882 

117.83376 

117.82825 

117.82402 

117.81625 

117.81154 

117.80243 

117.79672 

117.79195 

117.78524 

117.77847 

117.77298 

117.76727 

117.76056 

117.75384 

117.74819 

117.74353 

117.73799 

117.73083 

117.72261 

117.71690 

117.71 108 

117.70505 

117.69860 

I'T. 69322 

17.68801 
117.68041 
117.67486 

1  7.66764 

117.66293 

117.66094 

117.65706 

117.65518 

117.65447 

17.  65464 

117.65476 

117.65466 

j 32*65295 
117.65390 


32 

32 
32 
32 
32 
32 
32 
32 


48793 
48789 
49203 
50.:  49 
49304 
48779 
48750 
48722 
32.48723 
32.48722 
32.48717 
32.43740 
32.490  20 
32.48071 
32.48762 
32.48753 
32.48753 
32.48750 
32.48746 
32.4S741 
32.48683 
32.48717 
32.48709 
32.48697 
48587 
48680 
4870  3 
48697 

3 


32 
32 
32 
32 
32 
32 
32 
32 


32 
32 
32 
32 

32.4864 
32.48637 
32.48  6  60 
48o54 
43677 
48646 
48640 
48654 
48o28 
48594 
48560 
32.48550 
32.48602 
32.48606 
32.48600 
32.48736 
32.48728 
32.48579 
32.48573 
32.48565 
32.48591 
32.48692 
32.4909  1 
32.49386 
32.49745 
32.50144 
32.50738 
32.51333 
32.51884 
52.52479 
32.52818 

itmu 


ERR 

.u0035 

.00074 

-.COO  11 

.00208 

-.00454 

-.00290 

-.00119 

.20090 

.00  023 

-.00093 

-.00093 

-»  00:21*6 

-.00033 

-.00255 

-.00  302 

-.00003 

-.00  248 

.00  091 

-.00216 

-.00121 

-.00171 

-.00505 

0 

-.00265 

-.0016S 

Q 

.00025 

-.00194 

-.00102 

— .00040 

-.00  172 

-.00159 

.00036 

.00  002 

--00023 

.00043 

-.00093 

-.00006 

-.00041 

-.-J0009 

-,uC233 

.00079 

.00  153 

.00126 

.00136 

.00  158 

.0000  5 

-.000  85 

-.oooio 

.00064 
.00113 
.00131  - 

.U0007 
-.00  130 
.00192 
.00108 
.01052 
.00178 
.00133 
, 0  0  0  3  3 

.00050 


XX 

17.98735 
17.98533 
17.96151 
17.97800 
17.97644 
17.97056 
17.964  10 
17.95730 
17.94973 
17.94196 
17.95^62 
17.92591 
17.91778 
17.91--72 
17.9054  9 
17.896  20 
17.83914 
17.83262 
17.37o29 
17.869  74 
17.362  56 
17.35475 
17.84333 
17.84254 
17.83682 
17.83170 
17.82552 
17.81993 
17.81293 
17.80645 
17.600  65 
17.79449 
17.738  08 
17.78204 
17.776  13 
17.76990 
17.76343 
17.75756 
17.75183 
17.74629 
17.7401 1 
17.7331 1 
.17.726  6  3 
17.72041 
17.71426 
17.70799 
17.70203 
17.69645 
17.69003 
17.63396 
17.67743 
17.67165 
17.66735 
.17.66324 
17.66001 
17.65779 
17.65653 
17.65533 
17.65536 
7. 6547*3 
7  .  6'5444 


3  2 
32 

32 

32, 

32. 


YY 

32.48795 
52.48792 

32.43953 
4939^: 
49353 
4  9126 
43976 
4  3374 
32.43314 
32.48777 
32.43753 
52.48748 
32.48857 
32.43362 
32.48322 
32.43797 
32.48779 
32.43763 
32.43759 
32.48752 
32.43724 
32.48721 
32.4371 6 
32.43709 
32.43700 
32.48692 
32.48696 
32.43697 
32.43675 
3  2.4  3660 
32.4566' 
32.48658 
32.43665 
32.48658 
32. 4365r 
32.43644 
32.48633 
32.4362^ 
32.43596 
32.48577 
32.43587 
32.48595 
32.48597 
32.48653 
32.48683 
52.48642 
32.4  8614 
32.43594 
52.48595 
32.48630 
32.48816 
32.4904k 

32.49324 
32.49652 
32.500  87 
32.50535 
32.51105 
32.51654 

32.d2598 
32.53033 


SMOERR 

.000  35 
-.00307 
-.00344 
C  0835 
-.01218 
- .  0 13  61 
-.01292 
-.01154 
-.01216 
-.01359 
-.01334 
-.01633 
-.01543 
-.01494 
-.01426 
-.01442. 
-.01420 
-.012  17 
-.01276 
-.01269 
-.01 350 
-.01386 
-.01212 
-.01401 
-.012  83 
-.010  90 
-.01030 
-.00874 
-.00  675 
-.00347 
-.00010 
.00293 
.00622 
.00758 

•  00370 
.01114 
.01237 
.01201 
.01107 
.01 121 
.01096 
.01230 
.01221 
.01189 

•  0 1  3  5  5 
.01416 
.01333 
.01177 
.01197 
.010  97 
.01  19  3 
.01053 
.00902 
.006  3  9 
.00659 
.00521 
.00733 
.00353 

•  00 1 1 1 

.00290 
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VEHICLE  TWO 


117.74514 
1 17.74614 
117.74983 
117.75484 
117.75936 
117.75779 
117.75902 
117.76103 
117.76065 
117.76143 
117.76372 
117.76418 
117.76674 
117.76848 
1 17.76714 
117.77153 
117.77251 
1 17.77651 
117.77490 
117.77702 
117.77820 
117.77281 
117.77982 
1 17.78439 
117.78818 
117.79102 
117.79331 
117.79722 
117.79387 
117.79516 
117.80135 
117.80091 
117.80250 
117.80393 
117.80494 
117.81006 
117.81318 
117.81615 
117.81875 
117.82088 
117.82.417 
117.82278 
117.82312 
1 >7. 82.41  9 
117.82897 
117.83082 
117.85395 
117.85757 
117.83691 
117.83715 
117.83625 
117.83754 
117.84139 
117.84307 
1  17.84057 
117.84308 
117.83485 
117.84093 
117.84132 
117.84461 
117.84535 


Y 

32.5 

32.5 
32.5 


32 

32.5 

52.5 

32.5 

32.5 

32  »5 

32.5 

32.5 

32.5 

52.5 

32.5 

32.5 

32.5 

32.5 

32.5 

32.5 

32.5 

32.5 

32.5 

32.5 

32.5 

32.5 

52.5 

32.5 

52.5 

32.5 

32.5 

32.5 

32.5 

52.5 

32.5 

32.5 

32.5 

32  .5 

32.5 

32.5 

32.5 

32.5 

32.5 

32.5 

32.5 

32.5 

32.5 

32.5 

32.5 

5t  •  3 

32.5 

32.5 

32.5 

32.5 

32  .  5 

32  . 5 

32.5 

32.5 

32.5 

32  •  5 

32.5 

32.5 


8  60 
571 
050 

751 
879 
457 
3  1  0 

447 
451 
454 
448 
87  2 
723 
:  0  3 

726 

148 

996 

700 

566 

557 

698 

716 

4-j  9 

6  p  6 

672 

530 

574 

6  5  0 

53  1 

?75 

497 

615 

40  9 

463 

6  8  9 

33  6 

872 

523 

516 

50  7 

3  58 

6  0  9 

763 

0  42 

690 

843 

2360 

2   5 

2343 

22  23 

2510 

2619 

2931 

5351 

3531 

3422 

3650 

3671 

5623 

33  17 

4297 


17. 

"7 

I    . 

7. 

7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 

7 

•  . 

7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 

-7 

f  • 

7. 
7. 

7. 
7. 

"7 
t    « 

7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 

7. 
7 

7 

•  . 

7. 

7. 
7. 


7. 

-v 

I  * 
7. 
7. 
7. 
7. 
7. 
7. 
7. 
7. 
7'. 
7. 
7. 
7. 
3. 


)0G 
01000 
020 

0  3 ' 

30' 
05000 
06  00  i 
07000 
5  3  DOO 
39000 
10000 
11000 
120 
13000 
140  0  ) 
1500=* 

17000 
18000 
19  )00 

200  " 

2^■' 

22( 
23000 
240C  i 

2500  I 
260 
-7  000 
23000 
29  00... 
30000 
31  000 
-. 

33000 
3400- 
350 
36000 
37000 
5800'' 
3900C 

4 1 ", 

42000 

43000 

44000 

45000 

46000 

47(i 

4800 

49000 

500 

5 100, j 

52^  - 

53000 

54000 

-'  C.  ■  \  ~'\  -] 

5  c  000 
57000 

590  00 
00000 


X  X 

17.74514 
17.745  54 
17.74726 
17.750  29 
17.75392 
17.7  55  47 
17.75689 
17.75855 
17.7  59  59 
17.760  21 
17.76161 
17.762  64 
17.76428 

17.765  96 

17.766  4  4 
17.76847 
17.77009 
17.77266 
17.77355 
17.77494 
17.77624 
17.77437 
17.77685 
17.77987 
17.78319 
17.786  33 
17.789  12 

17.792  36 
17.79296 

17.793  84 
17.79685 
T7.79847 
17.80012 
17.30164 
17.80296 
17.80580 
17.308  75 
17.81171 
17.814  52 
17.81707 
17.81991 
17.82106 
17.82188 
17.822  81 
17.825  27 
17.32749 
17.83:i07 
17,83307 
17.334  6  1 
17.83562 
17.83587 
17.33654 
17. 838 4 S 
17.840  32  32 
1 f.  84042 


32  •  5 
3  2 . 5 
3  2.5 
32.5 

32.5 

32.5 

?  n      c 

— *  C-  •    _' 

32.5 
32.5 
32.5 
32.5 
32.5 
32.5 
52.  5 
32. 5 
32.5 

3  Z  .  3 

32.5 
32.3 
32.5 
3  2.5 
32.5 
32.3 
32.5 
32.5 
52.5 
32.5 
32.5 
32.5 
32.5 
32.5 
32 .  d 

5<C.3 

32.5 
32.5 
32.5 

52.3 
32.5 
5  2 . 3 
32.5 
3  2.5 
32.  5 
32.5 
52.5 
32.5 
52.5 
52. r 
52.3 
3  2.5 
32.5 
32.3 
32.3 

7.  C        C 


17 
17 
17 
17 
17 


84148 
8  3385 
83967 
84033 

3  4204 


17.84337 


32.5 

32.5 
32.5 
52.5 
32.3 
32.5 


1860 

1744 
18  67 
182C 

1844 
163  9 
1537 
1501 

1481 
1470 
1462 
1626 
1664 

177r> 
1522 
1711 
170  7 
1651 
1613 
1647 
1675 
1568 
1616 
1658 

153  5 
1470 

154  2 
1478 
143  6 
1461 
1522 
1477 
1471 
I338 
1469 
1630 
15  37 

;r59 

1538 
1653 
1638 
1638 

j  77  :■ 

2  >,2  5 
2017 

2140 
2179 
2  3  1  1 
2454 
2633 
2920 
3164 
3267 
3420 
3521 
5561 
3664 
3917 
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APPENDIX  3 
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VEHICLE  OUTPUT  DATA 


VEHICLE  ONE 


VEHICLE  TWO 


DIST 


117. 
117. 
1  17. 
117. 
117. 
117. 
117. 
117. 
117. 
1  17. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
1  17* 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 
117. 


9874 
9823 
9758 
9727 
9741 
9617 
9544 
9471 
9384 
9303 
9236 
9128 
9056 
9001 
8926 
8853 
8786 
8728 
8668 
8599 
8518 
8430 
8388 
8338 
8282 
8240 
8163 
8115 
8024 
7967 
7920 
7852 
7785 
7730 
7673 
7606 
7538 
7482 
7435 
7380 
7308 
7226 
7169 
7111 
7050 
6986 
6932 
6  880 
6804 
6749 
6676 
6629 
6609 
6571 
6552 
6545 
6546 
6548 
6547 
6539 
6539 


32.4879 
32.4879 
32.4921 
3  2.5005 
32.4930 
32.4878 
32.4875 
32.4872 
32.4872 
32.4872 
32.4872 
32.4874 
32.4902 
32.4887 
32.4876 
32.4876 
32.4875 
3  2.4875 
32.4875 
32.4874 
32.4868 
32.4872 
32.4871 
32.4870 
32.4869 
32.4868 
32.4870 
32.4070 
32.4864 
32.4864 
32.4866 
32.4865 
32.4868 
32.4865 
32.4864 
32.4863 
32.4863 
32.4859 
32.4856 
32.4855 
32.4860 
32.4861 
32.4860 
32.4874 
32.4873 
32.4858 
32.4857 
32.4856 
32.4859 
32.4869 
32.4909 
32.4939 
32.4975 
32.5014 
3  2.507  4 
3  2.5133 
32.5188 
32.5248 
32.5282 
32.5332 
32.5381 


17.0000 
17.0100 
17.0  200 
17.0300 
17.0400 
17.0  500 
17.0600 
17.0700 
17.0800 
17.0900 
17.1000 
1  7 . 1 1 00 
17.1200 
17.1300 
17.1400 
17.1500 
17.1600 
17.1700 
17.1800 
17.1900 
17.2000 
17.2100 
17.2200 
17.2300 
17.2400 
17.2  500 
17.2600 
17.2700 

17.2  800 
17.2900 
17.3000 
17.3100 
17.3200 
17.3300 
17.3400 
17.3500 

17.3  600 
17.3700 
17.3800 
17.3900 

17.4  000 
17.4100 
17.4200 
17.4  300 
17.4  400 
17.4  500 
17.4  600 
17.4700 
17.4800 

17.4  900 
17.5000 
17.5100 

17.5  200 
17.5300 
17.5400 
17.5  500 
17.5600 
17.5700 
17.5800 
17.5  900 
18.0000 


.2441 

117, 

.7451 

32.5186 

.2378 

117, 
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APPENDIX  C-l 
Exponential  Smoothing 

Exponential  Smoothing  is  a  special  kind  of  moving  average  that  does 
not  require  keeping  a  long  record  of  previous  positions.   Like  other 
moving  averages,  it  has  a  stable  response  to  changes,  but  the  rate  of 
response  can  be  adjusted  readily. 

Suppose  we  adopted  some  such  rule  as  follows:   To  get  a  new  estimate 
of  the  average  position  of  any  vehicle,  add  to  the  previous  estimate  a 
fraction  of  the  amount  by  which  new  positions  exceeds  that  estimate.   Note 
that  under  this  rule,  positions  below  the  estimate  is  the  same  as  saying 
it  exceeded  the  estimate  by  a  negative  quantity,  and  adding  a  fraction  of 
a  negative  quantity  would  decrease  the  estimate. 

The  fraction  used  is  called  a  "smoothing  constant"  and  is  represen- 
ted by  o^  .      We  could  abbreviate  the  rule  in  the  form  of  an  equation: 

New  estimate  =  old  estimate  +  o<  (new  position  -  old  estimate) 

or 

New  estimate  =    o(  (new  position)    +    (!-©()    (old   estimate). 
This   above  process   could  continue   indefinitely.      A   simpler   representation 
may  be   achieved   if  we  call   the  new  position  Position      (P   ),    the  next  pre- 
vious position  Position      (P-,),    the  position  two    sequences   ago   Position    , 
and   so  on,    with  the  position  k  sequences   ago  being  represented   by  Posi- 
tion,    (P   ).      Our  expression  now  becomes 

2  3 

New  estimate  =   <*  PQ   +    <*  (1   -  <*  )P1   +   <*  (1   -  cK  )    P 2   +    c*(l   -  oi  )   P 

k  k 

+    ...    +   o<(l   -  C(  )    pk+    (1   -  oi  )      (estimate  made  k 

sequences   ago) . 
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The  new  estimate  is  then  a  sum  of  positions  in  past  sequences,  each  multi- 
plied by  a  definite  weight.   If  the  sum  of  these  weights  is  equal  to  one, 
then  the  estimate  can  be  called  an  average. 

Therefore  the  basic  rule  of  exponential  smoothing  is 
New  average  =  <=*(new  position)  +  (1  -  o< )  (old  average). 

How  to  Select  c* 

The  value  chosen  for  the  smoothing  constant  «X  determines  how  much 
of  the  past  positions  has  any  significant  effect  on  the  estimate  of  the 
average . 

If  a  small  value,  say  o(  -   0.01,  is  chosen  for  the  smoothing  constant, 
the  response  will  be  slow  and  gradual,  since  it  is  based  on  the  average  of 
very  many  past  positions.   A  high  value,  say  o(   =  0.5,  will  cause  the  es- 
timates to  respond  quickly,  not  only  to  real  changes,  but  also  to  the 
random  fluctuations. 

To  evaluate  alternative  choices  of  a  value  for  the  smoothing  constant, 
the  accuracy  of  the  positions  must  be  considered  along  with  the  effects 
due  to  Auto -correlation.  \^J 
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Appi&roix  c-iii 

METHODOLOGY  TO  BE  USED  FOR  COORDINATE  CONVERSION  C47 


The  procedure  to  be  used  for  conversions  of  raw  Omega  data  to  abso- 
lute Latitude/ Longitude  (or  the  inverse)  will  be  based  on  linear  inter- 
polation between  a  set  of  appropriately  selected  and  stored  points. 
The  accuracy  of  the  procedure  will  be  insured  by  the  methodology 
used  to  select  the  points  to  be  stored. 

The  basis  for  the  conversion  in  either  direction  will  be  three  matrices 
.  of  stored  data.     The  first  two  define  the  crossing  points  in  latitude/ 
longitude  of  integral  (but  not  necessarily  consecutive)  lanes  of  Omega 
phase  differences.     These  points  can  be  defined  as: 

Ltj    =    longitude  of  intersection  of  Omega  lane  i  with 
Omega  lane  j 

£jj    =    latitude  of  intersection  of  Omega  lane  i  with 
Omega  lane  j 

where  the  first  index  "i"  refers  to  the  lane  number  assoc- 
iated  with  phase  difference  between  one  set  of  two 
Omega  signals  and  the  second  index  "j"  to  another. 

In  addition  to  the  position  in  latitude  /longitude,   it  is  necessary  to 
store  the  actual  Omega  lane  numbers  of  these  crossing  points* 


ftkr     =    lane  number  associated  with  the  ,rr"       stored  cross- 

>  th 

ing  point  for  the  "k"       station  pair.     There  will 

normally  be  only  two  station  pairs  being  considered 

at  one  time  and  it  will  be  assumed  that  Cl*r  >Q,k  r-x 

to  make  searching  more  convenient.     Note  that  for 
k=l,    r  =  i  and  for  k=2,    r=j. 

1.     Omega  to  Latitude  /Longitude 

Reference  to  figure  1  indicates  that  one  problem  to  be  solved 
is  to  determine  the  Latitude/ Longitude  (L*/£#)  of  the  Omega 
fix  (flf,    Q%)  where  the  subscripts  1  and  2  define  the  station 
pair  as  the  "k"  index  does  in  the  Q^r  matrix. 

By  some  search  procedure  it  can  be  determined  that  ft*  lies 
between  0^  „   and  ft1)B+1,    and  similarly,    Q3j  n  <  C%  <  a3/B+1. 


# 


* 


In  addition  Clx  is  closest  to  Qj^    and  Qa  is  closest  to  fls  £.     Th 


ie 


equations  defining  L«  and  I*  then  become: 
L*    =    LA*    +    ^•■^m+1  n  "  h&*&l.  An* 


1       (L»;n+1    -    L„(I1)     . 


Af2a 
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FIGURE  1:    Geometry  of  conversion  problem 
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